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Prehistory of lattice QCD 

Lattice Quantum Chromodynamics (QCD) is in its 41st year since the seminal paper of Kenneth Wilson in 1974, 
Confinement of Quarks [1].' From that time remarkable events took place in our understanding of strong interactions. 
The lattice definition of QCD is an important milestone in the development of particle physics theories. Thanks to the 
asympototic freedom and the reflection positivity lattice QCD is the QCD as opposed to the continuum formulation 
which is not well-defined mathematically beyond perturbation theory. 

Historically, it was the quark model that paved the way for the modern theory of the strong interactions. In 1963 
Gell-Mann and Zweig proposed that the low lying hadron masses may be described as many body representations 
of the SU{2) and SU{3) flavor groups. Later, Greenberg, Han and Nambu introduced the color degree of freedom 
postulating that each quark flavor transforms under the fundamental representation of the SU (3) color group and that 
all hadrons are color singlets. The deep inelastic lepton-nucleon scattering experiments at SLAC that took place in the 
end of 1960s revealed that quarks or partons, as Feynman called them, behave like free particles at high energies. It was 
then clear that any quantum field theory of strong interactions should have this property. Although non-Abelian gauge 
theories were known since 1954 their renormalizibility was shown only at the beginning of 1970s. This development 
led Gross and Wilczek as well as Politzer in 1973 to show that non-Abelian gauge theories are asymtotically free. 
Quarks are however confined into hadrons and this property was not shown until Wilson proved it for the strongly 
coupled lattice gauge theories. 


The lattice 

The Wilson formulation allows a fully non-perturbative treatment of lattice gauge theories. Yet, the expectation 
value of a typical observable of the theory is not known analytically. Lattice QCD in not an exception, eg. the three- 
dimensional Ising model and many more models in statistical mechanics and field theory have no analytical solutions. 
The only possibility that remains is to use Markov Chain Monte Carlo algorithms like the local heat bath algorithm 
that was devised by Creutz in 1979 [4]. In that calculation he was able to show that the string tension of the SU{2) 
and SU{3) lattice theories scales as expected by the asymptotic freedom demonstrating thus in principle that lattice 
gauge theories can be investigated by means of Monte Carlo simulations. A series of proposals for the calculation of 
the glueball as well as hadron masses followed (see Rebbi’s book for more details [5]). 


For a historical review of the lattice theory see the paper of Creutz [2]. For a more recent review see the one written by Ukawa [3]. 
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FIGURE 1. The running coupling constant of the SU{3) theory together with the one-loop and two-loop results (left plot) [7] 
and the light hadron spectrum of QCD (right plot) [9]. 


Moore’s law 

The first calculation of Creutz required around 1 MFlops to compute the plaquette of the SU{3) theory on a 6^ 
lattice, which could run on the VAX-11/780. However, the long distance effects of QCD require 30^ lattices at a 0.1 
fm lattice spacing, which means a GFlops machine. The quark-antiquark potential at those distances was calculated 
12 years later by Bali and Schilling on the CM-2 machine [6], whereas the running coupling of the SU{3) theory 
was calculated one year later by the Alpha Collaboration [7]. The calculation of the hadron spectrum without quark 
loops took 11 GFlops-years to 300 GFlops-years to complete on the GFll machine in 1993 and on the CP-PACS 
machine in 1999 respectively [8]. The full QCD spectrum was calculated in 2008 on the Blue Gene/P machine with 
the performance of around 100 TFlops [9]. Some results of these calculations are shown in Figure 1. It has taken about 
30 years from the first SU (3) calculation till the full QCD one. It is Moore’s law, the exponential increase of computing 
power with time, that enabled lattice QCD simulations to reach this milestone. 


Numerical algorithms 

QCD has not been solved yet. It will take some time until precise calculations of a large number of strongly 
interacting particle properties become pervasive. However, based on the great progress witnessed in the development of 
numerical algorithms, it is not expected to take another 30 years for the lattice QCD to achieve this goal. Looking back 
in time, the Interdisciplinary Project Center for Supercomputing (IPS) of the ETH Zurich was the ideal environment for 
a lattice group to develop numerical algorithms. Led by de Forcrand and influenced by a new version of the BiCGStab 
algorithm of Gutknecht, the Zurich group together with the Wuppertal group, led by Schilling in collaboration with 
Frommer, made a pioneering work on the non-hermitian solvers for lattice QCD, a research which resulted in large 
savings of computer time for quark propagator calculations [10]. From that time numerical algorithms for lattice 
QCD have attracted a pleiad of mathematicians including Gutknecht, Frommer, Golub, Saad, van der Vorst, Higham, 
Manteuffel and many more in a series of workshops of Numerical Analysis and Lattice QCD [11]. Among the recent 
achievements of this collaboration is the adaptive multigrid algorithm for the inversion of the Wilson operator [12]. 
When these ideas appeared, Liischer had already developed the concept of local coherence of the low lying modes of 
the Dirac operator resulting in a deflation based algorithm with large savings in computer time [13]. Figure 2 shows 
clearly that these algorithms accelerate inversions by a very large factor. 


Simulation algorithms 

Simulation of QCD with dynamical quarks has always been challenging. The Dirac operator must be inverted a large 
number of times along a trajectory generated by the Hybrid Monte Carlo algorithm [14]. The amount of computer 
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FIGURE 2. The local coherence of low lying modes of the Dirac operator is the basic property that allows an efficient deflation of 
the Dirac operator (left plot) [13]. Essentially the same algorithm may be constructed using the ideas of adaptive multigrid solvers 
(right plot, Babich et. al.) [12]. Both plots show a gain in computer time by more than an order of magnitude. 


resources using this algorithm at the physical point is of the order of Petaflop-years. The situation remained this 
way until Liischer made two important improvements; using a domain decomposition technique he could accelerate 
inversions by a factor of hve (see SAP+GCR performance in Figure 2) which translates to an accelerated Hybrid 
Monte Carlo algorithm by the same order of magnitude [15]. Deflation of the low modes of the Dirac operator [13] 
allowed him to devise the deflation accelerated simulation algorithm which gives an independent improvement by a 
factor five or so [16]. While the exponential increase in computing power for one invested monetary unit remains a 
very important driving force in lattice QCD, the algorithmic advances have been at least as important, particulary so 
in the last 15 years [17]. 


Chiral symmetry and the fifth dimension 

Chiral symmetry is an important symmetry of strong interactions. However, on the lattice, a Weyl fermion may not 
exist without its opposite chirality partner. This means that a Dirac fermion may not exist without its doubler. Wilson 
gave the doublers a mass proportional to the lattice cutoff destroying thus the normality of the Wilson-Dirac operator; 


[DwM^O. 


This way, one has to fine tune the bare quark mass in order to restore the chiral symmetry. It is easy to see that any 
75 -hermitian lattice Dirac operator is normal if it satisfies the Ginsparg-Wilson relation [18]; 

YsD + Dys = a DysD , 

where a is the lattice spacing. This simple relation has apparently only one formal solution, the Overlap Dirac operator 
of Neuberger [19]; 

aD = I + Dw{D*„Dwy^^^. 

Prior to this discovery a five dimensional theory of chiral fermions, the Domain Wall fermions, was formulated [20]. 
Since the Overlap operator is the Dirac operator of the four dimensional effective theory of Domain Wall fermions 
at any finite lattice spacing [21], lattice QCD may be regarded as an effective theory of a higher dimensional theory. 
Five dimensional theories are even more expensive to simulate and the hope is that various multigrid ideas will work 
in this case as well [22]. Another approach is to simulate local theories with the minimum number of doublers at the 
expense of a broken hypercubic symmetry [23]. Since the restoration of the broken symmetry involves the fine tuning 
of certain lattice operators [24] new versions of such fermions with manifest unitarity should be sought [25]. 









Balkan physicists contribution in lattice QCD 


Lattice physicists born in the Balkans and its close neighborhood include: I.O. Stamatescu (Heidelberg): RG flow. 
Topology, finite density and finite temperature. H. Neuberger (Rutgers): Chiral fermions, computation of the over¬ 
lap. Y. Shamir (Tel Aviv): Domain wall fermions, beyond SM. B. Svetitsky (Tel Aviv): Finite temperature, beyond 
SM. E.T. Tomboulis (UCLA): Topology, confinement. T. ^elik (Ankara): deconfinement, percolation. A. Vladikas 
(Rome): Non-perturbative renormalisation, heavy-meson decays. R Dimopoulos (Rome): Spectroscopy, pseudoscalar 
decay constants. C. Alexandrou (Nicosia): Heavy-light quarks, B-mesons, light baryons. H. Panagopoulos (Nicosia): 
Topology, finite temperature, perturbation theory. G. Koutsou (Nicosia): Spectroscopy, form factors. M. Constantinou 
(Nicosia): Form factors, non-perturbative renormalization. A. Tsapalis (Athens): Form factors, multi-quark potentials. 
K. Anagnostopoulos (Athens): Supersymmetric QM, large N super Yang-Mills. K. Farakos (Athens): Elektroweak 
pahse, gravity. P. Vranas (Livermore): temperature, domain wall fermions. K. Orginos (Williamsburg): Spec¬ 

troscopy, improvements, algorithms. D. Becirevic (Orsay): Spectroscopy, form factors, BB mixings. L. Levkova (Utah): 
Spectroscopy, decay constants, finite temperature. A. Alexandru (Washington): Spectroscopy, finite density, topology. 

S. Prelovsek (Ljubljana): Spectroscopy, scattering. M. Unsal (Raleigh): Topology, sypersymmetry, QCD-like theories. 

1. Hip (Zagreb): Ginsparg-Wilson fermions, topology. M. Marinkovic (CERN): Quark masses, topology, finite den¬ 
sity. D. Xhako (Tirana): Domain wall fermion algorithms. D. Djukanovic (Mainz): Spectroscopy, chiral perturbation 
theory, QED. R. Zeqirllari (Tirana): Minimally doubled fermions. T. Sulejmanpasic (Raleigh): Topology. L. Leskovec 
(Ljubljana) : Spectroscopy. M. Vidmar (Ljubljana): Spectroscopy. S. Zafeiropoulos (Frankfurt): Dirac spectrum. N. 
Bozovic (Wuppertal): Inversion algorithms. M. Blazenka (Zagreb): Decay constants. K. Hadjiyiannakou (Nicosia): 
Spectroscopy. A. Bori 9 i (Tirana): Numerical algorithms, Overlap/Domain wall fermions, minimally doubled fermions. 

The author would like to thank M. Creutz, Ph. de Forcrand and M. Liischer for their invaluable comments on this 
contribution. He is indebted to M. Marinkovic, M. Miiller-Preusker, B. Svetitsky and A. Vladikas for sharing their 
information on lattice physicists from the Balkans. 
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